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A STUDY OF MAN'S PHYSICAL CAPABILITIES 
ON THE MOON 

ABSTRACT 

A study was made to compare man's energy expenditure and gait 
characteristics, during self locomotion at various rates, in earth 
and in simulated lunar gravity conditions. Hie tests were made 
for the subject walking and running on the level and on grades up 
to 30° while in shirt sleeves and while wearing a suit pressurized 
to 3.5 psig. The results, presented in four volumes, may be useful 
for the design of space suits and life support systems and the 
planning of lunar exploration missions and their logistics. 
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FOREWORD 

This volume covers the investigation of lunar gravity simulation techniques in ”A 
Study of Man's Physical Capabilities on the Moon", NASA Contract NAS 1-4449. The 
prime contract was administered by Langley Research Center, with Mr. William Letko 
serving as NASA Technical Monitor. 

The investigation was performed by Case Institute of Technology under sub- 
contract to Northrop Space Laboratories where Dr. Walter Kuehriegger served as the 
Principal Investigator. Professor James B. Reswick, Director of the Case Institute 
engineering design center, guided and directed the investigation. The Cleveland 
Institute of Art provided extensive assistance to the Case Institute by fabricating the 
body harness and counterbalance suspension system, as directed by Professor 
Roy P. Hess. 

In view of the complexity and scope of the work performed under this contract, the 
final report has been organized in four separate volumes (numbered I thru IV). Since 
the work itself was broken down into phases it was possible to treat each phase individ- 
ually and document them correspondingly. The four volumes which comprise this report 
are identified as follows. 

CR-66115 Volume I, Part 1 - Lunar Gravity Simulation Facility 

CR-66116 Volume I, Part 2 - Instrumentation 

CR-66117 Volume II - Part 1 - Biomechanics Research Program 

CR-66118 Volume II - Part 2 - Biomechanics Research Program Appendices 

CR-66119 Volume III - Work Physiology Research Program 

CR-66120 Volume IV - Investigation of Lunar Gravity Simulation Techniques 

Volumes I thru III were produced by NSL and have been assigned Northrop Space 
Laboratories' document number NSL 65-153. Volume IV was prepared from material 
contributed by Case Institute and reports on their portion of the contract effort. The 
total report (all four volumes) summarizes the performance during the contract period 
from 2 November 1964 to 30 September 1966. 
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VOLUME IV - INVESTIGATION OF LUNAR GRAVITY SIMULATION TECHNIQUES 

By A. J. Bartley, R. P. Hess, R. J. Kirschensteiner, 

D. A. Millett, R. J. Morgen and J. B. Reswick 

Case Institute of Technology 
Engineering Design Center 


SUMMARY 


The opening sections of the report are concerned with the inclined plane method of 
simulating lunar gravity and some attempt has been made to improve the versatility of 
this method by refining the harness design. 

The bulk of the report deals with the development of a new method of simulating 
lunar gravity in which it is proposed to suspend the test subject in marionette-like fash- 
ion from a steel ceiling. A new concept of a magnetic air pad was used as the support at 
roof level and this enabled the subject to move freely with minimal frictional resistance. 
Negator springs were used to nullify five-sixths of the subject's weight. 

• The practicality of the system has been proved by building a 1/6 scale working 
model. A full-size system has not yet been built but sufficient research and design data 
has been collected to enable this to be done. 


INTRODUCTION 


In a previous report to Northrop the authors put forward a collection of new ideas 
in the field of lunar gravity simulation. The subject matter ranged over a wide field being 
concerned with body suspension systems, limb assist mechanisms, lunar topography 
simulation and improvements to the inclined plane type of simulator. In addition, a 
number of ideas were put forward with a view to developing some form of "vertical 
suspension" simulator which would be complementary to existing inclined plane techniques 
and which would enable the test subject to perform maneuvers which are difficult or im- 
possible to execute on the inclined plane. 

The final report is concerned with the practical exploitation of some of the more 
promising ideas in the previous report. The new work can be broadly divided into three 
sections: 


(1) Improvements for the inclined plane simulator 

(2) Design and development of a vertical suspension system 

(3) Design and development of body harness and limb support mechanisms 

In Section (1), the emphasis has been placed on the construction of a suspension 
system and harness which would give the test subject more freedom of movement than 
has been possible hitherto. 

In Section (2) consideration is given to the problem of developing a system in which 
the test subject would be suspended vertically from the ceiling of a room and yet be free 
to move in any direction within the confines of the room. The proposed solution involves 
the use of negator springs to give the desired uplift to the test subject and a magnetic-air 
type ceiling attachment to give the required mobile support at roof level. 

Finally, in Section (3) the design of harnesses and limb assist devices is con- 
sidered. This work has been done in connection with the vertical suspension system 
based on the magnetic air-pad and negator spring device. 
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SECTION 1 


IMPROVEMENTS FOR THE INCLINED PLANE SIMULATOR 


A number of improvements were applied to the Langley type inclined plane simu- 
lator (Ref 1) and these are illustrated in Figure 1. 

Lateral Movement of the Legs 

The support devices for the legs were attached to the same cable which went from 
each leg support up and over a pulley. This pulley was attached to the spreader bar 
which held the cables to all the supports for the body segments. 

This arrangement allowed the subject to spread and pull together his legs thereby 
achieving a more natural standing, walking, etc. , situation. Because of the mechanical 
arrangement of the cable and pulley the legs had to spread and pull together in unison; 
however, this duplicates normal leg action and was found to be a desirable feature, even 
though there was no need for a wide stance to maintain balance when landing from a 
jump, etc. 

Rotation of the Head 


A cable passing from a point on the front of the helmet at approximately the mid- 
point of the forehead and thence up and over a pulley to a point in the middle of the back 
of the helmet allowed azimuthal rotation of the head. No important actual area of free- 
dom was introduced with this device; however, it did cause the subject to feel less con- 
tained by the entire apparatus and consequently was considered to be a plus feature. 

Upper Arm Freedom of Swing 


A relatively large "C" brace was fabricated and its lower horizontal extension was 
attached to the upper body harness. The cable support to the spreader bar was attached 
to the upper horizontal extension of this "C M brace. This "C" brace replaced the cable 
which went directly from the upper body harness to the spreader bar, and in so doing 
eliminated cable interference when the upper arm was placed along side the body in a 
hanging down position or when the arm was swung backwards. Again no important func- 
tional freedom was introduced but the removal of another constraint was considered to 
be a plus factor by various subjects who tried the system with and without this innova- 
tion. 

Notes on the Development of an Additional Important Freedom of 
Movement for the Inclined Plane Simulator 

The freedom of rotation around the subject’s Z axis (Figure 2) seemed to be an 
important one to incorporate for it would allow the subject to turn round and re-trace 
his steps as and when he liked. 

As a result, thought was given to this problem and the resulting best solution in- 
dicated that two additional "C M braces extending in front of the subject would be required. 
One "C” brace would attach to the hip and chest supports and the other would attach to 
the present sling supported leg. 
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It was decided not to modify the existing equipment to include these additional "C" 
braces because there are a number of disadvantages associated with their presence. 
Among these are: 

1. The mechanism and bulk of two additional "C" braces in front 

of the subject appeared as if it would be psychologically inhibiting 
with regard to a sense of freedom on the part of the subject. 

2. The "C" braces extending in front of the subject would interfere 
with many crawling or bending movements. 

However, in order to obtain maximum versatility, it may be desirable to design future 
inclined plane simulators in a manner that would permit easy addition and removal of 
these extra braces. 

Inclined Plane System - Conclusions 

By suitable development and refinement, it is possible to build an inclined plane 
simulator which would give the test subject four + degrees of freedom and also provide 
him with a reasonably large test area in which to work. The experiments carried out so 
far have shown that it is possible to crawl, walk, run and jump with relative ease. In 
view of the simplicity of the apparatus the results were quite promising and there would 
appear to be little doubt that this type of simulator merits further development. 

Nevertheless, the inclined plane methods possess certain limitations. One of 
these is the lack of realism felt by the test subject who cannot completely ignore the fact 
that he is attempting to walk on a near-vertical wall instead of on the floor when using a 
short cable radius. This disadvantage can no doubt be partially overcome by suitable 
topographical design and psychological conditioning of the subject. 

.Another disadvantage is the difficulty of achieving movement along the subject's 
X axis which is almost vertical relative to earth. Assuming, as suggested previously, 
that it would be possible with relative ease to give the subject azimuthal rotation of at 
least 180° it would still be difficult to provide freedom of movement in the X direction 
because his weight would need to be counterbalanced with great precision to nullify the 
effects of one g earth gravity. 

For these and other reasons it was considered desirable to examine the practical 
possibilities of vertical suspension systems, and the work presented in the next section 
will provide some idea of the progress made in this field. 
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SECTION 2 


DESIGN AND DEVELOPMENT OF VERTICAL SUSPENSION SYSTEM 


In any type of lunar gravity simulator there are two major problems to be solved - 
namely, how to reduce the weight of the subject, and how to provide him with the maxi- 
mum amount of mobility. These two problems are interdependent to the extent that a 
solution to the one will usually limit the number of possible solutions to the other. 

In the vertical simulator now discussed, the aim was to suspend the test subject 
from the ceiling of a large room in such a way that he was able to move freely in any 
direction he chose by using a suspension system which incorporated a negator spring 
device and a magnetic air pad. Figure 3 gives some idea of the kind of system envis- 
aged. 

A Negator Spring System 


There are two types of negator spring units, both employing the same basic 
negator spring coil. These are the simple extension spring and the constant torque 
motor. 

The force obtained by unwinding a negator coil is almost perfectly constant (i.e. , 
force - deflection gradient is zero). Through various methods of manufacturing control, 
gradients ranging from slightly positive to slightly negative can be achieved. 

A dynamic analysis of the proposed system (see appendix - A) indicates that a 
negator having a slightly positive force - deflection gradient is desirable. The reason 
for this is that as the negator coil is unwound, an increasing force must be exerted to 
overcome the weight of the extended portion. 

The practical work on negator spring units has been confined so far to the design 
and manufacture of a 5/6 g negator unit for the prototype model shown in Figure 4. This 
consists of two simple coil springs mounted back-to-back to provide an uplift of 20 lb. 
over a displacement range of 30 inches. 

Generally it is considered that negator spring devices are well suited to lunar 
gravity simulator applications because of the many advantages they offer. These 
include : 

1. Constant force for large deflections. 

2. Low weight and compactness. 

3. Relatively low inertial effects compared with counter balance weight systems. 

4. Relatively little friction since there is no friction inherent in the spring. 

Magnetic Air Pad Development 

The magnetic air pad was an overhead support fixture. This device gave the sup- 
ported object three degrees of freedom. The first two degrees of freedom were hori- 
zontal translation in a plane parallel to the ceiling. The third degree of freedom was 


VOL IV 


7 




These various movements were ob 


rotation, either clockwise or counterclockwise 
tained with almost negligible friction. 


The device was structurally simple. It consisted of a flat disc with an air orifice 
at its center. Around the orifice a number of permanent magnets were arranged. An 
air line was attached to the orifice and a simple hanging fixture was incorporated in the 
disc for the attachment of the negator spring assembly. 



The magnetic air pad was placed in contact with an overhead steel plate and com- 
pressed air was forced through the central orifice. This caused the air pad to float at 
some small distance from the ferro -magnetic surface, thus giving the device its iow 
friction characteristic. 


A prototype device (see Figure 4) four inches in diameter fitted with an Alnico 
magnet and operating at 75 - 80 psi was found to support a load of 55 pounds demon- 
strating the feasibility of the method. 

Figure 5 shows a schematic of the magnetic air pad and gives the basic physical 
parameters which governed its performance. For static equilibrium in the vertical 
direction the following relationship was established. 


r F z = - f l + f m + 2 


R 

"f 


(P - P ) rdr 
v a r 


= 0 


or 


F L = F M + 27r 


R 

/ 


(P - P ) rdr 
v a r 


F t = F-, + F . 
L M air 


It is known that the magnetic force F^ is simply a function of the distance h. For 

any particular permanent magnet, this function can be determined by a simple 
experiment. 

The pressure acting on the upper side of the disc P^ was assumed to be a function 

of the radial location r, the height of the disc from the surface h, the stagnation pres- 
sure P q and the orifice diameter. Figure 6 shows the relationship between P r , r and h 

for fixed values of P and the orifice diameter. This data was obtained with the experi- 
o 

mental test arrangement shown in Figure 7. 

Theoretical works by Mori (Ref. 2) and Robinson (Ref. 3) explain the shape of the 
experimentally observed pressure profiles. These analyses predict a supersonic flow 
issuing from the orifice followed by a shock wave in the bearing clearance space. The 
analyses yield good qualitative agreement with the experimental results shown in 
Figure 6. 

The relationship between F . , the force due to the air pressure, and h, the air 

B.ir 

gap, can be determined directly from experiment or can be derived from the pressure 
profile data. Figure 8 shows the typical form of the relationship. The curve is divided 
into four regions - A, B, C and D. In region A the force is negative which means that 
the net force due to air pressure is in the negative z direction. In regions B and C the 
net force due to the air pressure is positive. In other words, an air pad without mag- 
nets is capable of supporting a hanging load in both regions. 
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(Ferromagnetic ceiling) 


XxXXXXXXXXXXXXXXXXXXXX 



r = Radius 

h = Air Gap 

R = Radius of the disc 

P^ = Stagnation pressure 

P q = Ambient pressure 


Figure 5. - Schematic of magnetic air pad 
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Pressure 






Figure 7. - Experimental test arrangement 

The slope of the curve in regions A and B was positive. This meant that the air 
pad operating in either one of these regions was stable.* For instance, if an operating 
point in region B was chosen, the air pad was in equilibrium at that point. If for some 
reason the air pad was perturbed such that h increased slightly, the force F^ in- 
creased, thus driving the air pad back to its original position. The same general argu- 
ment can be used to show that region A was stable. 

* 

Although these comments l’efer to static stability no problems have been encountered 
with dynamic instability and the chosen design has proved to have very good dynamic 
characteristics when compared to other configurations. See Gross, Ref. 4. 
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In region C and a portion of region D the slope of the F . versus h curve was 

air 

negative. Here, if the air pad was perturbed slightly such that h increased there was 
no driving force toward the original equilibrium position. These regions of negative 
slope were therefore unstable. 

The general shape of the magnetic force F M versus air gap curve is shown in 

Figure 9. It is worthwhile to note that the slope of this curve was negative. This 
meant that a magnet supporting a hanging load at some distance from a ferro -magnetic 
surface is by itself unstable. 

Recalling the derived equation for the static load which the magnetic air pas was 
capable of supporting, it is seen that the F^ versus h curve was simply the sum of the 

F . versus h curve and the F-. versus h curve, 
air M 

Figure 10 shows this summation. For stable operation, the air pad must operate 
on the positive slope portion, ABC, of the F^ curve. For fail-safe operation (in case 

of air supply failure), however, the pad must operate with an air gap less than h^ as 

defined by the intersection of the F . curve with the zero load point. In other words, 

air 

if h is less than h^, the effect of cutting off the air supply will be to cause the disc 

(plus negator spring assembly) to adhere to the steel ceiling by virtue of the magnetic 
attraction which is then greater than the load. 
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Air gap, h 

Figure 9. - Relationship of magnetic force vs. air gap 

Although the experimental and theoretical investigation of the physical behavior 
of the magnetic air pad is not yet completed, sufficient data is now available to permit 
the design of a system to support the weight of a fully grown man. 

Some experiments have been conducted on an annular orifice type pad and a 
central recess type pad but both designs proved to be inferior to the central orifice 
design employed in this application. 

Safety Devices 


Since the magnetic air pad is not rigidly attached to the ceiling some form of 
safety device will be needed in order to protect the user in the event of a break -away. 
This aspect of the problem has not been given much attention so far but is unlikely to 
produce any major difficulties. 

One feasible method is to suspend a catching grid immediately below the negator 
spring assembly. The grid could be supported by wires attached to the corners of the 
room at roof level or thereabouts so that when a break -away occurs the equipment is 
held in suspension well above floor level. 

The suspension wire system will need to incorporate spring loaded coiler-decoiler 
devices in order to take up and give out slack as the test subject moves within the con- 
fines of the test area. 

As soon as a break -away occurs, the coilers must lock — this can be achieved by 
fitting solenoid operated brakes to the coiler units. The sensing element which activates 
the brakes could be a microswitch attached to the magnetic air pad and making sliding 
contact with the ceiling in which case the switch is closed and the brake shoes are held 
in the open position. Once a break -away occurs, the microswitch opens and the power 
supply to the solenoids is shut off thus causing the brakes to lock on. 

The possibility of accidental operation of the microswitch by local variations in 
ceiling height can be obviated by suitable choice of microswitch travel — say, one-half 
inch. The system is basically "fail-safe" in operation since the brakes lock on when- 
ever the circuit opens irrespective of the cause. 
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SECTION 3 


DESIGN AND DEVELOPMENT OF BODY HARNESS AND 
LIMB SUPPORT MECHANISMS 


The body support system that has been designed applies 5/6 g negating forces to 
the body and limbs in a way which provides the most natural feel of 1/6 g simulation, 
along with the maximum possible freedom of body and limb movement and the minimum 
amount of interference with body and limb movement by supporting cables. In testing, 
it was found that this system has just under 4-1/2 degrees of freedom. 

The torso support system, which was prototyped, consisted of a rigid saddle 
structure which the subject straddled plus straps to and over the shoulders. His torso 
weight was distributed through four bone points on his pelvis; the anterior superior 
iliac spine below his waist and the ischial tuberosities below his torso. See Figure 11. 
The 5/6 g negation forces were taken through points at about the collar bone on the 
shoulder straps and at points on either side of the hips of the saddle structure. The 
harness was found to be very comfortable during all tests. Padding was provided at 
the shoulders and at the points of contact below the waist. 


Ischial 

tuberosity 


Anterior superior 
iac spine 



Figure 11.- Saddle structure 
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The inclined plane simulator was used as a control to compare the feel of 1/6 g 
simulation, freedom of body and limb movement and support interference. 

The lower leg support consisted of a brace running the length of the lower leg 
with an arm prejecting under the instep of the foot. This brace was fastened to the leg 
with straps above the calf and at the arch of the foot at the instep. The suspension 
cable attachment point was located below the pivot point of the knee. 

Support for the entire arm was provided by a brace strapped at a point below the 
elbow pivot point. 

The exact positions of the leg and arm suspension cable attachment points was 
chosen so that the weight of the upper limb was distributed equally between the body 
and lower limb supports. 

The negation load of the head was taken through a helmet which was fitted with a 
gimbal to allow forward, backward and side to side movement. The support cable 
allowed for pivoting. 

All support cables led directly to the negating units. Exceptions to this direct 
arrangement occurred in the leg support cables which went through guides attached to 
the body harness. These guides led the leg support cables to the back of the subject 
and out of the way of his arms for all controlled maneuvering positions. This resulted 
in some diminution of negation forces to the legs when the cable from the harness guide 
to the leg was not vertical but this was considered to be unimportant relative to the 
principal function of the guides. 

The present equipment was designed to give the body a pitch movement around 
the X axis of 110° and a roll movement about the Y axis of 15° in the vertical and prone 
positions. 

Experiments were conducted on the body and limb support system described above 
with the various segments of the body being counterbalanced by weights approximately 
5/6 of the weight of the segment. The supporting cables passed through ball bearing 
pulleys which were fastened to an overhead support 30 feet from the floor and then 
through another series of pulleys 12 feet from the first so that the hanging counter- 
weights would not interfere with the subject. The set of pulleys directly above the sub- 
ject was designed to pivot through 180° to enable the subject to move freely about the 
Z axis. The overall arrangement is sketched in Figure 12. 

Experiments on this vertical prototype mockup with casual subjects are shown in 
Figures 13 through 17. They indicated that the simulation of a 1/6 g situation was fully 
as good as on the inclined plane simulator. The vertical orientation of the subject 
heightened the feeling of reduced gravity because of orientation to normal surrounding 
objects. Gross-motor activities such as jumping, simulated walking, running or 
leaping had the same feel as on the inclined plane simulator. Falling over was almost 
as good as on the inclined plane. Support cable interference occurred when the subject 
rotated backwards around his X axis to simulate a rearwards fall. If the subject ro- 
tated backwards more than 15° from the vertical position, the leg and lower torso sup- 
port cables tended to pass around his arms to a position in front of him. The 15 
rearward rotation, however, was considered sufficient to indicate to the subject that 
he had lost his balance. 
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Support cable interference also occurred when the subject tried to make a com- 
plete crossover of his arms in front of his body. It is not known at present whether 
this interference can be easily overcome but if this semigross motor activity had to 
be executed, then the arm support cables could be disconnected. 

It was generally felt that the inertia of the counterweights caused movements to 
be somewhat slower than on the inclined plane, but this effect should be largely elimi- 
nated when negator springs are employed. 

Figures 13 through 17 are a series of photographs designed to illustrate the har- 
ness design and the freedom of movement available to the test subject. 
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CONCLUSIONS 


The work which has been carried out on the inclined plane simulator showed that 
this type of system could be further refined to give the subject additional freedom of 
movement and an attempt has been made to show how this could be done. 

The most important part of the report was devoted to the vertical suspension type 
simulator. This type of simulator will permit the test subject to perform maneuvers 
which are difficult or impossible to perform on the inclined plane. Although much work 
still remains to be done, the experience now available clearly shows that the method is 
practical. 

The system has the advantages that the test subject is supported in the normal 
upright position and is able to move anywhere within the confines of the test area. It 
also provides him with sufficient freedom to perform most normal activities such as 
walking, running, changing direction, etc. 

The present harness has been evolved after prolonged testing and experimenta- 
tion and now provides a high degree of comfort to the test subject. 

The next logical step in the overall development process is the design and con- 
struction of a full size prototype system; sufficient basic design data has now been ac- 
cumulated to enable this step to be taken. 
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RECOMMENDATIONS 


The following recommendations by tasks are made to indicate the importance of 
continuing this work in the interest of improving lunar gravity simulation techniques: 

Task 1: The work on the magnetic air pad systems should be completed by experi- 

mentation and analyses for recommended and optimal suspension units. 

These recommendations shall include the design, radius, the magnet and 
air bearing flow rates. 

Task 2: Should be performed to finalize the work on the negator spring system sus- 

pended from the magnetic air pad. This should include information for 
recommending the 

• selection of springs 

• attachment of springs 

• dynamics of spring including the fatigue life and replacement schedule 

Task 3: Design and development of a harness and limb support mechanism for verti- 

cal subject suspension. This harness should be made adjustable to accom- 
modate a range of subjects in the shirt sleeve and pressure suited condition. 

Task 4: Production of a vertical suspension prototype for the simulation of lunar 

gravity. This prototype should employ the recommended magnetic air pad 
developed under Task 1, the recommended negator spring system from 
Task 2 and the harness recommended in Task 3. A safety mechanism should 
be installed and tested frequently/ or prior to each test run. This safety 
mechanism must be designed to prevent at all times the magnetic air pads 
and remaining suspension gear from falling onto the subject or the floor 
below. 

Task 5: Tests should be performed with the prototype developed under Task 4. Two 

or more subjects should be used to determine 

• the physical space available for lunar gravity stimulation experiments 

• the types of experiments that can be conducted with the vertical suspen- 
sion device 

• the dynamic response limitations, distortions or any other limitations 
observed during the tests 

• the effects of partially negated suspension versus total suspension 
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APPENDIX A 


ANALYSIS OF SYSTEM DYNAMICS 


Introduction 


This appendix contains an abbreviated version of an analysis of the dynamic 
behavior of a vertical suspension type simulator The system to be analyzed is shown in 
Figure Al. and is a simplified version of the real system shown in Figure 3. 

The system components consist of one magnetic air pad, a back-to-back negator 
spring assembly and the subject mass. The subject mass will be considered to be inani- 
mate. The negator spring assembly housing is attached to the air pad by a ball joint. If 
the system is confined to move in a plane perpendicular to the ceiling, its position can be 
defined by three generalized coordinates; namely, the position of the air pad center of 
gravity, the angle between the extended negator springs and the vertical, and the distance 
from the ball joint to the subject mass. An important assumption must be made, however, 
before these three coordinates can correctly determine the position of the system. This 
is that the extended portion of the negator spring does not bend in the plane of motion. 

This is probably a valid assumption for the orientation of the negator springs shown in 
Figure Al.* If, however, the negator spring assembly were rotated through 90 degrees, 
the width of the negator springs would be perpendicular to the plane of motion. Due to the 
very low rigidity of the negator springs in this direction, the assumption would no longer 
apply and the three coordinates would not suffice in locating the system. 

The following are the assumptions which will be made in the derivation of the 
dynamic equations: 

1. The extended portion of the negator spring is rigid and inextensible. 

2. The air pad moves with negligible friction. 

3. The up-and-down dynamics of the air pad are so small that they can be 
neglected. 

4. The ball joint has negligible friction. 

5. The air drag on the system is negligible. 

6. The motion of the system is planar. 

7. The negator spring system has a constant force for any elongation. 

The parameters used in the analysis are defined as follows: 

x Location of the center of gravity of the air pad. 


* Note that the width of the negator spring is parallel to the plane of motion. 
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6 Angle that the extended portion of the negator spring makes with the vertical. 

I The distance from the ball joint to the center of gravity of the subject mass. 

a The distance from the center of gravity of the air pad to the ball joint. 

b The distance from the ball joint to the center of gravity of the negator spring 

spools and housing. 

r Radius of a negator spring spool. 

s 

R Radius of the coiled portion of negator spring. 

W Width of the negator spring, 
t Thickness of the negator spring. 

L Total length of a negator spring. 

D co . The center of gravity of the system of masses suspended from the ball joint 
& measured from that point. 

X A distance parallel to any point on the extended portion of the negator spring, 
x Horizontal distance from a reference point, 
y Vertical distance from a reference point. 

Mass of the subject. 

M Mass of the air pad. 

r 

M^ Mass of the negator spring housing. 

M Mass of a negator spring spool containing bearings. 

M Mass of the extended portion of one negator spring. 

V 

M Mass of the coiled portion of one negator spring. 

P The density of the negator spring material. 

I Moment of inertia of one spool about the ball joint. 

I Moment of inertia of the coiled portion of one negator spring about the ball joint. 

*ho Moment of inertia of the negator spring housing about the ball joint. 

I Moment of inertia of one spool about its axis of rotation. 

I Moment of inertia of the coiled portion of one negator spring about the axis of 

rotation. 
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F The force output of both negator springs. 

T The kinetic energy of the system. 

Q A generalized force. 

Damping coefficient. 
t A torque. 

Tf Coulamb bearing friction torque. 

Kinetic Energy of the Extended Portions of Negator Springs 


With reference to Figure Al, the position of any element of mass along the 
extended portion of the negator spring can be written 

x, = x + A sin 0 (Al) 

A p 

y^ = -a - A cos $ (A2) 

The component velocities of any element of mass can be found by differentiating 
the above positions with respect to time. 

x^ = Xp + X sin 0 + A 0 cos 6 (A3) 

y^ = -X cos 6 + A 0 sin $ (A 4) 

The kinetic energy of the extended portions of both negator springs can be written 
as 

l 

T e = Pwt / (x A 2 + y/) dA (A5) 

b 

Upon substitution of the squares of the component velocities into the above expres- 
sion, the kinetic energy can be written 

l 

T g = pwt j (x p 2 + X 2 + A 2 0 2 + 2x A sin $ + 2x p A ftcos^) dA (A6) 

b 

It should be noted that X is a constant over. the integration which takes place at any 
instant of time. In addition, X must be equal to l due to physical constraint. Upon com- 
pleting the integration and making the above-mentioned substitution, the kinetic energy 
of the extended portion of both negator springs can be expressed as follows: 

T g = pwt {[x p 2 + j 2 + 2x p i sin 0] [/- b] + [j^ - b 3 ] + x p 0cos 6[ ^ - b 2 ]} (A7) 
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Kinetic Energy of the Subject Mass 


For the sake of simplicity, the mass will be assumed concentrated at a point. 
The coordinates of this point are 


x T = x + l sin 0 
L p 

A^ = -a - l cos 6 


(A8) 

(A9) 


The velocities of this point are found by differentiating the above coordinates with 
respect to time. 

(A10) 


x T = x + /sin $ + i 6 cos 0 

Li p 

y L = - /cos0 + /Osin 9 


(All) 


The kinetic energy of the mass can be written 

t l = i M L (*L 2 + yL 2 ^ 


(A 12) 


Upon making the appropriate substitutions, the kinetic energy of the subject mass 


is 


T t = \ M t (x n 2 + / 2 + / 2 £ 2 + 2x /sin 6 + 2x / 0 cos $ ) 

Jr r 


L 2 L v p 
Kinetic Energy of the Negator Spring Housing 

The coordinates of the center of gravity of the negator spring housing are 

x, = x + b sin 6 
h p 

y, = -a - b cos 0 


(A13) 


The velocities are 


x. = x + b 9 cos 6 
h p 

y^ = b 0 sin 9 


(A 14) 
(A 15) 

(A 16) 
(A 17) 


The kinetic energy of the housing can then be written as 

Tu = (x 2 + b 2 6 2 + 2x b(9cos 0) + 

h 2 h ' p p 2 hcg 


(A18) 


Using the parallel axis theorem, the moment of inertia of the housing about the 
ball joint is 


l un = 1. + M. b 

h^ hcg h 


(A19) 
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Substituting this equation into the expression for kinetic energy yields the following: 


T h = 5 M h( i p 2 + V &0O8 ») + | 1 h/ 


(A20) 


Kinetic Energy of Both Negator Spring Spools 

The coordinates of the center of gravity of the two-spool combination are 


x = x + b sin 0 (A21) 

s p 

y = -a - b cos 9 (A22) 

s 

The corresponding velocities of this center of mass are 

x = x + b0cos0 (A23) 

s p ' 

y = b 9 sin 9 (A24) 

s 


The spools rotate with an angular velocity a . This angular velocity is a function 
of both l and l, and can be expressed as 


<*= | R - ^(L +b -£) + r g 2 (A 2 5) 

where R is the variable radius of a coil. The expression for R is derived in the next 
section. 

Solving for the kinetic energy due to translation of the center of mass plus that due 
to rotation about the center of mass gives 


T = M [x 2 + 2x b0cos 0] + I 9 2 + I a— - ^ 

s s L p p J so s r t 


P P 

Kinetic Energy of the Coiled Portions of Both Negator Springs 


[77 (L + b - ^ ) + r g ] 


The coordinates of the center of mass of the coil are 


x = x + b sin 9 
c p 


y = -a -b cos0 
c 

It should be noted that these coordinates are not a function 
of the variable l. 

The velocities of the center of mass are 

x = x + b 9 cos 0 
c p 


(A 26) 


(A27) 

(A28) 


(A29) 


32 


VOL IV 



I 


y c = b(9 sin 0 


(A30) 


The kinetic energy of the coils can be expressed as 


• 9 

T c * M c (i p + 2i p b « cos e ) + 1 0# e 2 + i oc 

R 


The parameters M , I a and I are functions of the distance /. 

C CC7 etc * 

It can be easily seen that the mass of one coil is 


M - pwt (L + b - l) 


(A31) 


(A32) 


The radius of the coil can be found from the evaluation of the following integral. 

R C 

/ rdr = f ~27 T dS ( A33 > 


length. 


Where C is the length of the coiled negator spring and ds is an element of arc 


C = L +b -l 

Upon evaluation of the above integral, R is found. 


R= v/£(L+b-/)+r s 2 

The moment of inertia of one coil about its cylindrical axis is 


I pw 

ac 2 


„4 4 

R ~ r s 


(A34) 


(A35) 


(A36) 


The moment of inertial of a coil about an axis perpendicular to the axis of 
symmetry through the ball joint is 


I = m , 
6c c L 


< r2 - r s 2 L + W S 


+ b 


■1 


found. 


4 12 

Making the appropriate substitutions, the kinetic energy of both coils is 


(A37) 


T = pwt (L + b - £ ) (x 2 + 2x b 9 cos $ ) 
c p p 
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+ p wt (L + b - l) 



(L +b -l 



.2 

e 


[±r(L+b-£) + 2r 2 ] i 2 
+ | pwt (L + b - l) 

[-^(L+b-/) +r g 2 ] 


(A38) 


Kinetic Energy of the Air Pad 

The kinetic energy of the air pad can be simply expressed by 


ml - 2 
T =tt m x 

p 2 p p 


(A39) 


Kinetic Energy of the System 

The kinetic energy of the system is the summation of the kinetic energies of 
the component masses. 


T = t + T„ + T, + T +T +T 
e l h s c p 


(A40) 


Generalized Forces 

Because of the negligible air drag assumption, the generalized force due to a 

variation in x, is 
b 



Using the assumptions of negligible air drag and friction in the ball joint, the 
generalized force due to the variation of 0 is 

Q 0 = " g (M L + M h + 2M g + 2M c + 2M e ) D cg sin 6 

Where D is the center of mass of the system suspended below the ball joint, 
eg 


(A41) 


D = 
eg 


b (2M g + M h ) + /M l + 2pwtbL + pwt(i - b)' 

M t + M, + 2M + 2M + 2JV1 
L h s c e 


(A42) 


can now be expressed as 


- - g [ b (2M + M T ) + iM T + 2pwtbL +pwt(f - b) ] sin0 (A43) 
s Ij 
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In calculating the generalized force due to a variation in St , the frictional torque in 
the bearings of one negator spring spool is assumed to be of the form 

T = - (D a + -SL T) (A44) 

m | a | 1 

The frictional torque contains a viscous friction term and a coulomb friction term. 

It will be assumed that these are independent of the load applied to the bearing because 
the variation in the bearing load should be small. 

2r 

.The frictional force due to this torque can be written f - ~ ~ making the substitution 
• A At 


The frictional force due to the bearings in the spools is 


1 I ~i 
i = -2 (D + -pr 

R 2 \l 


R 


)• 


The generalized force due to a variation in l is 


Q i = (2M e + ) S cos 6 ~ 2 (D m 


R 


M 


T f> 


(A45) 


(A46) 


Substituting the expressions for and R into the above equation, is found in 
terms of the generalized coordinates. 

2 [D m i+ TiiW* <L + b ' /, + r s 2] 

Q . - [ 2 p wt (/ - b) + M t ] g cos 0 - t s F (A47) 

1 L (L + b -. / ) + r g ^ 


Derivation of Dynamic Equations 

Using the method of Lagrange, three dynamic equations will be derived. There will 
be one equation for each generalized coordinate or degree of freedom. 


The first equation of motion for the system is 


Then ~- 
dt 



) = 0 


d , _aT_ _ er = 0 
dt 1 ax ’ ax xp 


dr 

ax 


= o = 



(A48) 
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It is seen that the coordinate can be ignored from the standpoint of the dynamics 
of the system. Furthermore, the above equation is immediately integrable and has the 


dT 

0 * 


const. 


Since the quantity — — is the momentum associated with the coordinate x , the 

P 

previous equation is a statement of conservation of momentum in the direction. 
Carrying out the required differentiation, this equation becomes 

(Mp i + 2M g + 2 pwtL + M^) Xp + (M^ + 2M g + 2 p wtL) b 0 cos 0 

+ ( 1 sin $ +£ $ cos $) + pwt [ (/ - b) 2 0 cos 0 

+ 2 (/ - b) / sin 0] = const. (A49) 


The second equation of motion for the system is 

d dT _ JTT 

dt 00 3 0 

Performing the required differentiations yields: 


= Q 0 


(A50) 


0 


2/3 If 3 ~ b3 ) + V 2 + + 2I s# + ( l h h - l ) 2 


2 pwt (l b - -jg - + b 2 ) 


+ 61 


2 pwt/ 2 + 2 /M l 


2 1 r 

b -ij - 2 pwt ^ + b 2 j + Xp pwt + b 2 ) cos 0 


1 M^/ cos0 1 2M g b cos0 + 2pwt^L - /jb cos 0 


+■ it 


b^2M g +■ * /M l + 2 pwtbL + pwt (f - b) 2 


sin 0 0 


(A. - ) 1 1 
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The third equation of motion for the system is 


dt \ g] 

The evaluation of this equation yields: 


d / dT \ _ dji 
\ dJ / dl 


Q, 


1 


pwt (l- bj 


*■ M t + 21 + 

L s a 


pwt (l + b -i>jj -jj— (l + b - + 2r g 


+ r 


x^ | 2 pwt sin - bj + sin #J 

j_ i 2 I w ( L + b -<) 4~[( L + b -^) t2r 3 2 l l 

\w ( L+b -') + r s 2 f 
(L+b-/) + -yy— +b 2 ] -pwt< 2 -M L <| 


•2 r t 

+ e pwt 


jspwt (i - bj + j g cos 9 - F 


n + r 1 \ I 4 ( L ' b ~t )+ r 2 1 

D m + T f f7TV«^ l sJ 


^(L+b-l) 


+ r 


(A 52) 


(A53) 


The three equations of motion are highly nonlinear and contain a great degree of 
dynamic coupling. 


Linearization and Solution o f Dynamic Equations 

An attempt will be made to linearize and solve the previously derived equations of 
motion. Due to the highly nonlinear nature of the equations, it is obvious that any linear 
approximation will have an extremely limited application. 

The First Equation of Motion 

By making the following limitations and assumptions, the first equation of motion 
can be linearized: 


9 ~ is limited to small angles 
9 - is small 
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x - is small 
P 

P. - is small 

The generalized coordinate P can be written l # = X + € where € is small and A 
is constant. Differentiating i with respect to time l -- € . 

Substituting the above into the first equation of motion and using cos 6 = 1 and 
sin 6^6, the equation becomes: 

(M + M, + 2M + p wtL + M T ) x + (A54) 

' p h s H V p 

(M, + 2M +2 P wtL) b(9 M t [ k 0 + (X + € ) 0 1 + 
ns ij 

pwt [ (X + € -b)^ '6 + 2 (X+ € -b)€ 6} - const. 

Dropping terms of order two or more in the above-mentioned parameters, the 
first equation of motion is linearized. 

• • 

Ax - B 6 - const. 

P 

where A and B are the following constants: 

A = M + M, + 2M + 2 p wtL + M T (A55) 

p h s H L 

B - (M h + 2M g + 2 pwtL) b + M L X + pwt [ X - b ] 2 (A56) 

The Second Equation of Motion 

Assuming the following limitations: 

$ - is limited to small angles 
'0 - is small 
# - is small 

and dropping terms of second order or more, the second quation of motion takes the 
form: 

6>C + x D + 0E + 6 e F + x € G = 0 (A57) 

P P 

where C, D, E, F and G are constants. 

If it is now required that the accelerations 0 anc * * x * are small and again drop- 
ping terms of second order, one obtains a linear differential equation. 

0C+x D + 0 E - 0 
P 
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where 


C - 2/3wt(A 3 -b 2 ) + M t A 2 + I, + 21 . 

Li no s o 

2 2 

+ ^r£ (L + b-A) + 2pwt(L + b-A)(^ + b 2 ) (A 58) 

D = pwt(A 2 + b 2 ) + M L \ + M h b + 2M g b + 2 pWt(L - A)b (A 59) 

E = g[b(2M + M T ) + AM T + 2 p wtbL + wt (X - b) 2 ] (A60) 

S JU Ij 

Solution of the First and Second Linearized 
Equations of Motion 

Differentiating the first equation with respect to time: 



Substituting 

obtained 


into the second equation, the following equation in 


6 is 


<9 + [ 


C 



6 = 0 


(A61) 


From physical considerations, it is evident that 9 will be a periodic function, 
and therefore the solution of the above equation will be sinusoidal. 

Using the initial conditions 


0=0 at t = 0 
0 = a sin pt 


where 



The period of motion will be 


(A62) 


T 


2n 

P 
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Now that the motion of 6 is known, the motion of the air pad can be derived. 

,, W •• • • 

Remembering that x = - -r 6 and making the substitution for 0 

P A 

Xp = p 2 sin pt (A63) 

Using the following initial conditions 

x = v n at t = 0 

P 0 

x = 0 at t = 0 

P 


and integrating twice with respect to time, the following equation for x^ is obtained: 


Ba 

x - - -7— sin 

P A 


inpt+(v 0 +^)t 


(A64) 


It is observed that the frequency of oscillation for both x and 0 is the same. 

B P 

The amplitude of x , however, will be -r times greater than the amplitude of 6. 
p A 


It should be stressed that these derived functions are only as valid as the lineariza- 
tion assumptions. If they do have any physical relationship to the real system, it will be 
for very short periods of time and small motions. 


The Third Equation of Motion and Its Solution 

By imposing the following limitations and assumptions, the third equation of motion 
can be linearized. 

6 - is limited to small angles 

6 - is small 

x - is small 

P 

R - is constant 


The assumption that the outer radius of the coiled portion of the negator spring, R, 
is constant is a good one since the negator spring is very thin. As a direct consequence 
of this assumption, the moment of inertia of the coiled negator spring can be written 


I ca = M c R 2 = pwt (L + b -/ ) R 2 


(A6 5) 
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Making the above-mentioned assumptions, the third equation of motion has the 
following linearized form: 


i' +C 1 l +C 2 l + C 3 = 0 


(A66) 


where the constant coefficients have the following values: 


2D 


C 1 = 


m 


1 1 (2pwtL + M l ) + 2l ga 


C 2 = 


~(2pwtg) 


21 


2pwtL + M + 

L R 2 


sa 


[2pwtbg-M L g F +,4i ^ 


C 3 = 


171 


21 


2pwtL + M T + „ 

L R 2 


sa 


It should be noted that the coefficient changes value with a change in sign of l . 

This means C„ will have one value while the subject mass is rising and another during 
the fall. 6 

2 

When a and b are two real distinct roots of the characteristic equation S + Cj S 
+ C 2 = 0, the differential equation has the following solution: 


i(t) 



1 1 at 1 bt 

C 2 a (a - b) e b (b - a) e 


■\i 


«>V 


c, i 


( 0 + ) 


at 


a - b 


b - a 



+ / ( 0 + ) 


a 

a - b 



b 

b - a 


bt 


e 


(A67) 


The constant D , which is the damping coefficient, is found to be very small for 
m ^ 

roller bearings lubricated with light oil. If D m is sufficiently small, the constant 

can be neglected in the differential equation. For the case where is neglected, the 
following equation of motion is obtained: 


l(t) 



j cosh at + / (0 + ) ~ 
a 


sinh at 


(A68) 


where a = | y~C 2 | and is real. 
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Differentiating with respect to time, the velocity and acceleration are calculated. 

c 

i (t) = a |/(u + ) - — ^ | sinh at + J i (O') J cosh at (A69) 


a 


l (t) = a* 1 1 (0 ) - — | | cosh at + a | i (0 + ) | sinh at 

a 


(A70) 


The time at which i (t) = 0, which is the instant that Cg changes value, is 


t = 
P 


- tanh -1 
a 


a 


jjjOJ- 


i(0 ) - -y 


(A71) 


where C„ is evaluated at t . It can also be shown that the acceleration l will have a 
3 • p 

minimum at the time t . 

P 

For the purpose of judging the performance of this lunar gravity simulator, l (t) 
should be investigated. For good simulation, this acceleration should approximate l/6g 
or the gravity on the moon. 

Sample Calculations 

In the following calculations, the following parametric values will be assumed: 


M 

P 

= 

1. 56 x 10 ^ slugs 

m l 

= 

. 558 slugs 

M h 

= 

3 x l(f 3 slugs 

M 

s 

= 

_3 

4. 5 x 10 slugs 

2pwt 

= 

5. 8 x 10 4 slugs/in. 

P 

= 

. 00879 slugs/in 3 

w 

= 

1. 5 in. 

t 

= 

. 022 in. 

R 

= 

1.75 in. 

L 

= 

144 in. 

A 

= 

120 in. 

b 

— 

1. 80 in. 
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J sa 

= .0208 slug-in? 


-3 

= 9.7x10 slugs-: 

-3 

1 se 

= 14.6x10 slugs 

D m 

= 3. 2 x 10 ^ in. -lb- 

T f 

= .010 in./lb 

g 

= 32. 2 ft/ sec 2 


2 

n. 

. 2 


Calculation of Up-and-Down Motion 


For the purpose of this calculation, a subject mass of . 558 slugs was chosen. 

This mass will have an earth weight of 18 pounds. As a first attempt at designing a 
simulator, negator spring force of 15 pounds is chosen. With this spring force, almost 
5/6 of the subject mass’s earth weight will be negated statically. 

The subject mass is now given an initial velocity upward and its motion is observed. 
The initial conditions are 


then 


[ (0+) = 120 in. 
i (0 + ) = -88 in. /sec 

C 2 = -• 342 sec -2 

2 

Cg = -54. 1 in. /sec 
a = -. 5845 sec 1 


Due to the fact that the damping coefficient D is very small, can be 
neglected and the following equation applies 


£(t) = -7T + 
C 2 


+ ^3 

no ) - -4 

a 2 J 


• _|_ 

cosh at + 1 (0 ) — sinh at 
' ' a 


The rise in time is found to be t = 1.032 sec. At this time the subject mass will 
have zero velocity and a distance 1 of 76. 3 in. , which corresponds to a lump height of 
43. 7 in. . For the fall, the value of Cg will change due to the influence that the change 

in direction of velocity has upon the friction in the bearing. The new constant is 

2 

Cg = -54. 5 in. /sec . The drop time is then found to be t = 1. 032 sec. 

In general, the drop and rise times will not be the same. However, under this 
set of conditions, there is negligible difference between them. 
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The height above the floor and the acceleration of the subject mass for this jump 
are shown plotted versus time in Figures A2 and A3 respectively. For the purpose of 
comparison, the values of height and acceleration for a jump on the moon under the 
same initial conditions are also plotted. 

It is evident from these graphs that the simulation is far from perfect. Namely, 
the jump height on the simulator is 43. 7 in. compared to 59.9 in. on the moon. In 
addition, the acceleration is still much higher than the lunar acceleration. It can be 
concluded from this that the choice of a 15-pound negator spring which will give good 
static simulation is not adequate dynamically. 

To demonstrate this point, a new negator spring force of 15. 8 lb is chosen. 

The height and acceleration versus time curves for this spring force are also plot- 
ted on Figures A2 and A3. As can be seen, the simulation is greatly improved, 
but as would be expected, still not perfect. It is important to note that the accel- 
eration has an almost hyperbolic shape rather than being constant. 

Calculation of Oscillatory Motion 


Keeping in mind the large number of restrictions necessary for the solution of 
the first and second equations of motion, the following analysis is presented. 

Solving the equations derived in the section on the third equation of motion for 
0 and Xp about the operating point A = 120 in. , the following is obtained: 

The angular frequency p = 2. 12 rad/sec 

The period T = 2.96 sec 

and the equations for 6 and Xp are 

6 = a sin (2. 12t) 

Xp ~ —31. 8a sin (2. 12t) + (67.4a+V Q )t 

Recalling that the total jump time for the 15. 8 lb negator spring was 2.64 sec- 
onds, it is evident that, if the above equations have any physical significance, it will 
be only for a very small portion of the total period. 

Conclusion 

The equations of motion of a lunar gravity simulator were derived using the 
method of Lagrange. These equations were found to be highly nonlinear; and, 
therefore, an attempt was made to linearize them. Some of the assumptions neces- 
sary for the linearization are of a questionable nature. This fact must be kept in mind 
when reviewing the results of these equations. 

It was found that, for the up-and-down motion of the subject mass, a negator 
spring force which would give good static lunar simulation would not yield good dynamic 
characteristics. It was also found that by adjusting the negator spring force to a higher 
level, the dynamic characteristics could be improved. 
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Figure A2. - Comparison of height vs. time curves for an initial velocity of 88 in. /sec 
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Figure A3. - Comparison of acceleration vs. time curves for an initial velocity of 
88 in ./sec 
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Recalling that it is possible to obtain negator springs with a positive force gradi- 
ent, it is probable that such a spring would flatten out the acceleration versus time 
curve for a jump and therefore improve the dynamic simulation. 

The influence of a positive force gradient negator spring and the possibility of 
numerical solutions to the nonlinear differential equations will be the objects of 
future work. 
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